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The multiparticle entanglement in the Lipkin-Meshkov-Glick model has been discussed exten- 
sively in this paper. Measured by the global entanglement and its generalization, our calculation 
shows that the multiparticle entanglement can faithfully detect quantum phase transitions. For an 
antiferromagnetic case the multiparticle entanglement reaches the maximum at the transition point, 
whereas for ferromagnetic coupling, two different behaviors of multiparticle entanglement can be 
identified, dependent on the anisotropic parameter in the coupling. 
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I. INTRODUCTION 



The exploration of the connection between statistical 
mechanics and quantum information has been extensive 
in recent years since the work [l| . Especially the research 
of entanglement in many-body systems has contributed 
to the comprehensive crossover between the two hot ar- 
eas 0. Furthermore, the finding of integer or fractional 
quantum Hall effect in two-dimensional many-body sys- 
tems imposes a challenge on the universal understanding 
of phase transitions, since the traditional theory for phase 
transition cannot incorporate these novel phenomena [|| . 
Recently the research of quantum entanglement in two- 
dimensional many-body systems provides the clear char- 
acterization for different quantum orders [3| ■ These facts 
suggest that quantum entanglement would play a vital 
role in the understanding of many-body effects. 

Bipartite entanglement was first studied, and fo- 
cused on the connection to the criticality in spin-chain 
systems [j| @. This interest comes from the fact that 
quantum phase transition is related to the construction of 
the long-range correlations in many-body systems. Hence 
it is a natural conjecture that quantum entanglement, 
as a depiction of the non-local correlation, could detect 
the appearance of long-range correlation. Great progress 
has been made for the block entanglement in many-body 
systems; the area law of block entanglement entropy has 
been generally constructed by the conformal field theory. 
Furthermore the violation of the area law has been iden- 
tified as a reliable detection of quantum phase transitions 
in one-dimensional systems(see Ref. [2] for a comprehen- 
sive review). However the situation becomes complex for 
high-dimensional systems: the violation of the area law 
in one-dimensional case when the system is critical, does 
not seem to hold in higher dimensions Even for pair- 
wise entanglement in many-body systems the results are 
not satisfying . For example, the cutoff in the defini- 
tion of concurrence may lead to unphysical results when 
one focuses on the connection of entanglement and phase 
transition in many-body systems Q. 
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The situation becomes more complex for multiparti- 
cle entanglement (ME) because of the absence of unified 
measurement for ME[8(. However, it is a natural specula- 
tion that ME should play a more fundamental role than a 
bipartite one for the understanding of many-body effects 
with consideration of the universal interaction in many- 
body systems. Recently the discussions of ME in many- 
body systems have been given more attention because 
of the availability of some special entangled states, e.g., 
cluster states for one-way quantum computation[9[, n- 
party Greenberger-Horne-Zeilinger(GHZ) state and W- 
state[l(|. Although great effort has been devoted to the 
measurement of ME, the analytical or operational mea- 
surements exist only for some special cases Q. The con- 
nection of quantum phase transition and ME has also 
been discussed extensively in [lf|, [H[ ■ However these 
discussions are mainly on the one-dimensional spin-1/2 
XY model, and the difficulty of calculating ME obstructs 
further exploration. 

Recently, the global entanglement has been con- 
structed for the quantification of ME by Meyer and 
Wallach[l4j], which possesses the virtues of the avail- 
ability of analytical expression and operability. More- 
over global entanglement is measurable experimentally 
since it is directly related to the mixedness of single 
party 15[ . Another important character of global en- 
tanglement is the monotonicity under local operations 
and classical communication (LOCC), if one notes that 
global entanglement is intimately related to the linear 
entropy [HJ. Consequently, Oliveira and his collabora- 
tors improved this definition for measuring some special 
entangled states, e.g. n-party W state or GHZ state|13|. 
Moreover the connection of the generalized global en- 
tanglement and quantum phase transition has also been 
explored in the one-dimensional spin-1/2 XY model, in 
which entanglement reaches the maximum near to the 
transition point fl3|. 

It is interesting to note that the nearest neighbor cou- 
pling is beneficial to the formation of ME in the one- 
dimensional spin-1/2 XY model. Since the particle cor- 
relation is short-ranged in this model one should note 
that the maximum of global entanglement maybe come 
from the distribution of pairwise entanglement [l2l [l7| . 
Hence it is tempting to present a discussion about ME 
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when the correlation is long range and the coupling is be- 
yond the nearest neighbor case. Fortunately the Lipkin- 
Meshkov-Glick(LMG) model 18] provides us the bench- 
mark for exploring this point since the collective interac- 
tion in this model. Then it is expected that ME would 
play a critical role. 



Recently the entanglement in the LMG model has been 
extensively studied, such as concurrence IT9I. [20L [2lT| . one- 
tangle [12], entanglement entropy (23l 12411251] and gener- 
alized entanglement [26]. The concurrence in the LMG 
model displays sensitivity to the appearance of quantum 
phase transition [ill [2l|], except for some special cases [2lj|. 
A possible explanation of this discrepancy is that the 
trace operation performed in the calculation of concur- 
rence inevitably kills some correlations between spins [2 lj. 
With respect to this point, Barnum and his collaborators 
constructed a subsystem-independent measure of entan- 
glement, based on a distinguished subspace of observables 
for the svstem[27j. The named generalized entanglement 
introduced by Barnum, et.al. has also been discussed in 
the LMG model, which displayed the ability of detecting 
the phase transitions ,26] . Moreover the authors show the 
equivalency between generalized entanglement and the 
global entanglement defined by Meyer and Wallachflij]. 
However, as shown in [26| . it is indispensable for the con- 
struction of the distinguished subspace of observables to 
obtain the knowledge of the ground state in many-body 
systems, that in most cases is very difficult. The research 
of entanglement entropy in the LMG model shows that 
the entropy was divergent under thermodynamic limit 
near the phase transition point, and moreover shows a 
discontinuity at the critical point for the isotropic cou- 
pling case (23. 



The generalized global entanglement (gGE), defined 
by Oliveira, et. oL[13|. is a generalization of the global 
entanglement (GE). With respect to the equivalence be- 
tween Barnum's generalized entanglement and GE, gGE 
provides a universal characterization of ME in many- 
body systems. Hence, it is interesting to present a com- 
prehensive research of gGE and GE in the LMG model. 
Our discussion also presents detailed research for anti- 
ferromagnetic coupling and some interesting results can 
be obtained, which is rarely touched on in the previous 
works. I should point out that the goal for this paper fo- 
cuses on the connection between ME, measured by gGE 
and GE respectively, and quantum phase transition in 
the LMG model. For this purpose, the paper is orga- 
nized as following. In Sec. II the Hamiltonian is pre- 
sented, and ground states are determined analytically. 
The phase diagram will be identified by introducing the 
proper parameter. In Sec. Ill the analytical expressions 
for gGE and GE are presented. Based on these formulas 
the multiparticle entanglements for ferro-magnetic and 
antiferro-magnetic couplings are discussed respectively. 
The conclusions and discussions are given in Sec. IV. 



II. HAMILTONIAN AND GROUND STATE 

The LMG model describes a set of spin-half particles 
coupled to all others with an interaction independent of 
the position and the nature of the elements. The Hamil- 
tonian can be written as 



H = -±-{Sl + 1 Sl)-h z S z 



N 



(1) 



in which S a = ^ i=1 a l a /2(a. — x,y,z) and the a a is 
the Pauli operator, and N is the total particle number 
in this system. The prefactor 1/N is essential to en- 
sure the convergence of the free energy per spin in the 
thermodynamic limit. Anti-ferromagnetic or ferromag- 
netic interaction can be obtained dependent on A < or 
not (A ^ 0). The Hamiltonian preserves the total spin 
and does not couple the state having spin pointing in the 
direction perpendicular to the field, i.e. 



JV 



[tf,S 2 ]=0, [H,Hai]=0. 



(2) 



For isotropic coupling 7=1, [H,S Z ] = and the spec- 
trum of Eq.([T]) can be determined exactly. However for 
7^1, the spectrum can be determined in principle by 
Bethe-type equations [28j and the analytical expressions 
are difficult to obtain. 

A distinguished character of Eq. ([T|) is the collective in- 
teraction, which is the same for any particle and indepen- 
dent of the space configuration of the system. Because of 
long-range correlation between particles, the mean-field 
analysis is adaptive for this model[29j]. The research of 
phase transition in the LMG model has shown that there 
is a second-order transition at h = h z /\\\ = 1 for the 
ferromagnetic case and a first-order one at h = for the 
antiferromagnetic case[l9l. [2^]. 

A proper parameter for characterizing the phase dia- 
gram is the total spin in the direction z for the ground 
state. For ferromagnetic coupling, one has 



1 - 2{S Z )/N = 



0, 



h > 1 



1-h, he [0,1), 



(3) 



which corresponds to the disorder-order transition, and 
obviously the point h = 1 is a second-order phase transi- 
tion point. This phase transition could be attributed to 
the disappearance of the energy gap; at the symmetric 
phase h > 1 the energy gap above the ground state is 
finite, whereas at the broken phase < h < 1 the en- 
ergy gap vanishes under thermodynamic limit [20]. For 
antiferromagnetic coupling, 



2(S Z )/N = 



1. 



h > 
-1, h< 0: 



(4) 



Obviously there is a first-order phase transition at the 
point h = 0. For this case the energy gap above ground 
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state vanishes only at the transition point h = un- 
der thermodynamic limit, and no level crossing happens 
when h £ o[2l|. 

The ground state for 7^1 can be determined analyt- 
ically with the help of Holstein-Primakoff(HP) transfor- 
mation and low-energy approximation [2jj. In Ref. [30j ] . 
the ground state has been obtained with the considera- 
tion of the finite number effect. The general expression 
reads 



Iff) 



, [N/2] 
c £ — ' 

n=Q 
[N/2] 



(2n- 1)!! 
2nH 



tanh"ar|2n) 



EH)' 



(2n- 1)!! 
2^!! 



tanh 2 " x 



(5) 



in which |2n) is the Fock state of the boson opera- 
tor introduced by Holstein-Primakoff transformation and 
[N/2] denotes the integer part not more than N/2. One 
should note that the determination of the ground state 
Eq. (O is based on HP transformation, which preserves 
the symmetry Eq. @, and the following discussion is 
heavily based on this ground state. Dependent on the 
style of interaction, tanh _x has different expressions. For 
ferromagnetic case A > 0, it satisfies the relation(l9j 
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tanh 2x = 



h>l 



2/1-1-7 ' 



(6) 



For antiferromagnetic coupling A < 0, it is determined 
by 



tanh J2x 



1 
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l + 7+2|/i|' 



(7) 



For isotropic case 7 = 1, the calculation is exact. The 
ground state can be formulated generally as \g) = \S = 
S z = M) . For ferromagnetic coupling, 



M 



I[hN/2], 0<h<l 
h > 1. 



(8) 



in which I[n] expresses the integer not more than n. For 
antiferromagnetic coupling , 



1-f, h<0. 



(9) 



intimate connection to bipartite entanglement [la | , it has 
been proven that the global entanglement is operational 
and more importantly, monotonic under LOCC. A sim- 
plified ex pre ssions for global entanglement is provided by 
Brennen [la ] 



N-l 



2 (! - jf E Tr ^l)- 



(10) 



k=0 



For the LMG model, one can obtain Q(\g)) = 2(1-Trp 2 ), 
in which pi stands for the single-particle reduced den- 
sity operator. Furthermore, Oliveira and his collabora- 
tors have improved this definition in order that it can 
measure some special entangle states, e.g. ® n \EPR) n or 
n-party GHZ state. The main procedure is to measure 
the entanglement between any two parties and the oth- 
ers, and then average all possible bipartition[l3j]. In the 
LMG model, for the symmetry of particle permutation, 
the generalized global entanglement can be written as 



E g = -{l-Tv[p 2 2 ]). 



(11) 



in which pi denotes the reduced density operator for any 
two particles. p\,pi can be determined through the cor- 
relation functions 13111 



(<Tl a cr 2 /3) 



4(Sg)-iV 
N(N-l) 
2{[S a ,Sp] + )-N 
N(N — 1) 



(12) 



in which a, (3 = x, y, z. 

With respect to the symmetry Eq. {2} and the ground 
state Eq. ([S]), one can obtain GE and gGE respectively 



E g = l-h(2(a i 



Q(\g)) = i-{<r 

> 2 + (<JlxV2x) 2 



\0\z02z 



f|3) 



Based on Eqs. © and ifTZ)). ME in the LMG model can 
be decided analytically, and some interesting properties 
can be found. The discussion below is divided into two 
cases: one focuses on the anisotropic coupling, and the 
other is for isotropic coupling for which the exact results 
can be obtained. 



III. MULTIPARTICLE ENTANGLEMENT 



Recently, Meyer and Wallach have constructed the 
global entanglement for measuring ME in spin systems. 
The main procedure is to first measure the entanglement 
between any party and the others, and then calculate the 
average of all possible bipartition[Ij|. Although the crit- 
icism that it is not a genuine ME measure because of the 



A. anisotropic coupling 

The analytical results can be obtained under large N 
limit with the hypothesis that the excitation would only 
happen for the low energy states[li|. Based on Eqs. © 
and p2p . one obtains 
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[JV/2] 



n— 

[iV/2] 



7V(iV- 1) 



tt o V - 2n + 2)(jV - 2n + l)2n(2n - l)c 2 „_ 2 c 2n + 2n(jV - 2n)cJ„] 

1 r z — ' 



-2 



71=0 

[JV/2] 



iV(iV- 1) 



p-2 J2 - 2n + 2)(iV - 2n + l)2n(2n - l)c 2 „_ 2 c 2n - 2n(jV - 2n)<4 



n=0 



[JV/2] 



N(N- 1) 



—2 ^ 2n(n - 2n)c^ 



(14) 




FIG. 1: (Color online) The multiparticle entanglement for 
ferromagnetic coupling, measured by gGE (denoted by E g 
and solid lines) and GE (denoted by Q(\g}) and dashed lines) 
vs the rescaled magnetic field h. We have chosen 7 = 0.5(a) 
and 7 = 0(b) for this plot. The green and black solid lines 
correspond to N — 50, 500 respectively. 



in which c 2 „ = (~1)" y ^ 2 2 „(p !! tanh" _x and tanhx is 
decided by Eqs.© and ©. From Eq. |(T3J), ME can be 
determined analytically. 

-Ferromagnetic case- gGE and GE have both been 
plotted for different 7 in Fig[TJ It is obvious that ME 
reaches the maximum closed to phase transition point 
h = 1 and the slope of curves tends to be infinite. In 
recent papersJUJ, the authors have shown that the sin- 
gularity of gGE is directly connected to the degeneracy 
of the ground-state energy at the phase transition point. 
Since the energy gap above ground state vanishes at crit- 




FIG. 2: (Color online) E g (black triangles) and Q(\g)) (green 
triangles) for ferromagnetic coupling vs. the particle number 
N at the phase transition point h = 1. One should note that 
tanh 2a; for any 7 £ [0, 1) is identical in this case. 



ical point h = 1, the singularity of GE and gGE can be 
attributed to the degeneracy of ground-state energy, and 
could be used as a reliable detector for the phase tran- 
sition in this case. Furthermore the finite-size scaling 
at the phase transition point h = 1 displays the non- 
sensitivity both of gGE and GE to the particle number 
N, as shown in Fig[2j 

-Anti] err omagnetic case- The situation is intricate. We 
have plotted GE and gGE vs. h by choosing 7 = (a) 
and 7=1/2 (b) respectively, in FigO Since there is 
a first-order quantum phase transition, the figures show 
that GE and gGE both are maximum at transition point 
h = 0. However a further calculation shows two different 
behaviors for ME; For 7 G (0,1) GE and gGE show a 
cusp at phase transition point h — 0, which means that 
the first derivation of gGE and GE with respect to h is 
discontinued but finite at h = 0. Since the degeneracy of 
ground-state energy happens only at h — 0, this discon- 
tinuity of gGE and GE is attributed to the degeneracy of 
ground-state energy, whereas for 7 = 0, the figure shows 
that gGE and GE both have a drastic increase closed to 
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FIG. 3: (Color online) E g (solid lines) and Q(\g}) (dashed 
lines) for antiferromagnetic coupling vs. the rescaled mag- 
netic field h. The parameter 7 = 0.5(a) and 7 = 0(b) have 
been chosen for this plot. Since the system is invariant under 
the changing of h «-> —h, these plottings are only for h > 
with N=50 (green lines) and N=200 (black lines) respectively. 



the phase transition point, and the first derivations of 
gGE and GE with h tend to be divergent. Furthermore 
our calculation shows that with the increasing of particle 
number, gGE and GE decrease for 7 £ (0, 1) at the tran- 
sition point, as shown in Fig. [U Moreover, GE and gGE 
have similar behaviors and the difference between them 
is slight. 



B. isotropic coupling 

When 7 = 1, the exact results can be obtained. With 
respect to Eqs. ([5]), © and (p~5|) . one has in this case 



_ 8M 2 

Q{\g)) = i-O 2 - 



2(4M 2 - Nf + {AM 2 - N 2 f 



QN 2 {N-lf 



(15) 



-Ferromagnetic case- With respect to Eqs. ([8]) and 
(fT5| . both gGE and GE are zero for fa > 1, indepen- 
dent of the particle number N since the ground state is 
a direct-product state of N particles with the same spin 
orientation from Eq. ([8]). As shown in Fig. [5j gGE and 
GE both are continued under thermodynamic limit at the 
phase transition fa = 1. This behavior is different from 
the conclusion made in Refs. 120,0,^3, in which entan- 
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FIG. 4: (Color online) E g and Q{\g)) for antiferromagnetic 
coupling vs. the particle number N at first-order phase tran- 
sition point h = 0. We have choosen three representative 
values of 7 for this plot. 
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FIG. 5: Eg (solid line) and Q(\g)) (dashed line) vs h for 
isotropic ferromagnetic coupling with N — > 00. 



glement entropy and concurrence both are discontinued 
at fa = 1 under thermodynamic limit. The main rea- 
son for this discrepancy is stated below. One notes that 
Eas. (fT5)l are a function of M/N. It is obvious from Eq. 
© that M/N is continuous at the transition point fa = 1 
under thermodynamic limit. Hence, it is not strange that 
gGE and GE are also continuous. Comparably the con- 
currence in [20] is redefined by adding a prefactor N — 1 
to keep finite under thermodynamic limit. Similarly for 
the calculation of entanglement entropy, all possible bi- 
partition has to be considered in the calculation to keep 
the entropy finite under thermodynamic limit (24j . which 
I point out plays the same function of the prefactor N — 1 
for the calculation of concurrence. Under the thermody- 
namic limit, the prefactor would play a nontrivial role. 
Since the calculation for entanglement has been imple- 
mented respectively in different regions because of Eq. 
([8]), the difference, which should disappear under ther- 
modynamic limit, may become finite because of this pref- 
actor. While, since our definition of ground state Eq.© 
has naturally considered the finite-number effect and GE 
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and gGE are the functions of the correlations, one does 
not need this prefactor to keep the measurements of en- 
tanglement finite under thermodynamic limit. Together 
with respect that the equivalency between generalized 
entanglement and GE have been proved [26j. the measure 
gGE may also has the great virtue of independence on 
the concept of subsystem. 

-Antiferromagnetic case- With respect to Eqs. §§§ and 
([T"5"|) . gGE and GE both vanish independently on N. 
Since the states for M = ±N/2 are degenerate at phase 
transition point h = 0, the ground state is undoubtedly 
the superposition of states \N/2,N/2) and \N/2,-N/2), 
written on the basis of {S 2 7 S z }, with equal weight, 

\g) = 4=(W2, N/2) + \N/2, -N/2)). (16) 

Then in this case, a genuine maximally multiparticle en- 
tangled state, so called n-party GHZ statefioj. can be ob- 
tained at the point ft = for a finite particle number, and 
under thermodynamic limit, it corresponds to the cele- 
brated ( Schrodinger cat ) macroscopic quantum super- 
position state. Obviously, the measurement of entangle- 
ment is discontinued at ft = 0, where a first-order phase 
transition happens under thermodynamic limit (32| . 

IV. DISCUSSIONS AND CONCLUSIONS 

Some comments and discussions should be presented. 
In this paper an extensive discussion of multiparticle en- 
tanglement, measured by GE [l4[ and gGE [13j], is pre- 
sented in the Lipkin-Meshkov-Glick model. Our discus- 
sion focuses on two different situations: for ferromag- 
netic coupling A > 0, when the anisotropic parameter 
< 7 < 1 gGE and GE both reach the maximum at the 
second-order phase transition point ft = 1, as shown in 
Figs. [1] Moreover they are nonsensitive to the variation 
of particle number N, as shown in Figf2] Whereas for the 
isotropic case 7=1 gGE and GE are zero at the phase 
transition point, shown in Fig. 

Another important situation is the appearance of an- 
tiferromagnetic coupling A < 0, for which there is a first- 
order phase transition at transition point ft — 0. gGE 
and GE both are calculated, as shown in Fig. [3] It is 
interesting that some different behaviors can be found in 
this case; one is that gGE and GE have a cusp at phase 
transition point when 7 € (0, 1), which means that the 
first derivation of gGE and GE with external magnetic 
field is discontinued but finite at h = 0, shown in Fig. 
[3Ja). Another case happens when 7 = 0, in which both 
gGE and GE have a drastic changing closed to ft- = 0, 
shown in Fig. [3](b) . Moreover our calculations show that 
gGE and GE for 7 € (0,1) decrease with the increment 
of N as shown in FigUJ whereas for 7 = they are non- 
sensitive to the particle number N . It is more interesting 
for the isotropic case that the entanglement is vanishing 
for ft ^ and has a sudden changing at h = 0, where 



a genuine maximally multiparticle entanglement state, 
n-party GHZ state for finite particle number, can be ob- 
tained. With connection of a scheme of the realization 
of the LMG model in optical cavity QED[33|, this result 
provides a powerful method to create ME experimentally. 

It was naturally expected that ME should be maxi- 
mum at the phase transition point since the correlation 
between particles would be long-range because of the ap- 
pearance of critical quantum fluctuation at the phase 
transition point. However an exceptional case appears 
in our discussion, which happens for ferromagnetic and 
isotropic coupling. In my own opinion, a reason for the 
difficulty in constructing the connection between entan- 
glement and quantum phase transition is that the up- 
to-date measurements for entanglement are generally a 
nonlinear function of correlation functions in many-body 
systems. Hence the singularity of correlation functions 
may be canceled [3J|. As a concrete illustration one notes 
that for antiferromagnetic coupling, M has a discontin- 
ued change at ft, = for 7 = 1, as shown in Eq.©. 
However from Eq. (|15|) it is obvious that the discontinu- 
ity of M for ft > and ft < has no effect on ME since 
gGE and GE are the functions of M 2 . 

Regardless of this defect, some interesting information 
can be obtained from the research of ME. An interest- 
ing speculation from our discussion is that the different 
finite-size scales may show the different state structures 
for the entanglement at the phase transition point. As 
shown previously in Ref. [131 ], with the increment of 
particle number the measures for entanglement for some 
states are decreasing, whereas for other states tend to 
be steady values. Since the increment of particle num- 
ber, or more generally the degree of freedom, means the 
stronger correlation between the particles in many-body 
systems, one can conclude that some entangled states are 
immune to the effect imposed by the increment of corre- 
lation between particles, whereas others are sensitive to 
the changing of correlation. With respect to the pursuit 
of decoherence-free space [35], our discussion may provide 
some useful information. 

With connection to the researches of the bipartite en- 
tanglement in many-body systems, one could note that it 
is difficult to construct a universal classification of phase 
transition based on the entanglement and its derivation. 
The main obstacle, in my own opinion, is the absence of 
physical definition of entanglement, i.e., how and what 
to define an " entanglement operator" . Since quantum 
entanglement is an important physical resource and can 
be measured experimentally, I believe in the existence of 
this operator. Fortunately a few works have attributed 
to this direction [26l. [36| . I hope our discussion will help 
in the understanding of entanglement in many-body sys- 
tems. 
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